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PREFACE 
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Arnold  Engineering  Development  Center  by  the  University  of 
Cincinnati,  Department  of  Aerospace  Engineering  under 
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were  B.  N.  Srivastava,  M.  J.  Werle,  and  R.  T.  Davis.  The 
Air  Force  Project  Engineer  for  this  contract  was  E.  R. 
Thompson,  AEDC/DYR.  The  Program  Element  Number  was  6580 7F. 

This  report  covers  the  work  during  the  period  April 
1974  to  March  1976.  The  reproducibles  used  in  the  repro¬ 
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I .  INTRODUCTION 


The  reentry  problem  continues  to  provide  motivation  for 
study  of  low  Reynolds  number  high  Mach  number  blunt  body  flows. 
This  flow  regime  is  one  in  which  the  viscous  effects  influence 
a  significant  portion  of  the  total  shock  layer  thickness, 
thereby  violating  the  classical  boundary  layer  approximations 
and  requiring  the  use  of  a  more  comprehensive  set  of  governing 
equations.  Such  a  set  are  the  viscous  shock  layer  equations 
(Refs.  1  and  2)  which  represent  an  intermediate  levef  of 
approximation  between  the  boundary  layer  and  Navier-Stokes 
equations.  These  equations  contain  all  of  the  terms  in  the 
Navier-Stokes  equations  which  contribute  to  second  order 
boundary  layer  theory  plus  those  which  arise  to  second  order 
in  the  outer  inviscid  portion  of  the  shock  layer.  Although 
some  evidence  exists  that  the  viscous  shock  layer  model  will 
suffice  for  reentry  type  flow  problems  (Refs.  3  and  4) ,  there 
has,  as  yet,  not  been  a  critical  assessment  of  the  range  of 
validity  of  these  equations.  This  situation  is  due  to  the 
difficulties  involved  in  solution  of  these  equations.  Although 
several  methods  have  been  presented  for  solving  the  "thin" 
shock  layer  approximation  to  these  more  general  equations 
(Refs.  5  and  6),  such  approaches  suffer  two  limitations.  First, 
they  are  based  on  the  assumption  that  the  pressure  gradients 
normal  to  the  body  surface  are  established  entirely  by 
centrifugal  effects,  and  second,  that  the  shock  wave 
lies  parallel  to  the  body  surface.  In  an  attempt  to  remove 
these  limitations,  methods  have  been  developed  (see  Refs.  7 
and  8)  for  addressing  the  full  shock  layer  equations  through 
a  relaxation  process  wherein  the  thin  shock  layer  assumptions 
are  removed  by  iteration.  However,  no  attempt  has  been  made, 
as  yet,  to  assess  the  range  of  validity  of  these  equations. 

This  is  performed  in  the  present  work  through  comparison  with 
experimental  stagnation  point  data  for  spherical  nose  shapes. 
Numerical  solutions  of  the  viscous  shock  layer  equations  are 
obtained  by  combining  an  implicit  finite  difference  scheme 
with  a  relaxation  technique  for  determining  the  bow  shock 
shape  (see  Ref.  7 ) •  The  effects  of  thin  layer  approximations, 
wall  slip  and  shock  slip  boundary  conditions  are  included  in 
the  analysis  to  establish  their  relative  importance.. 


II.  GOVERNING  EQUATIONS  FOR  THE  VISCOUS  SHOCK  LAYER 


The  governing  equations  are  written  in  a  boundary- layer 
coordinate  system  (see  Fig.  1) .  The  equations  and  notations 
used  are  the  same  as  those  used  by  Davis  (Ref.  1  )  or  Srivastava, 
Werle  and  Davis  (Ref.  9  ). 

Continuity  Equation; 

[(r  +  n  cos((i) '^pu]  g  +[(1  +  kn)  (r  +  n  cos<ti)'^  ~ 
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s-Momentum  Equation; 


p{u[Ug/(3.  +  kn)  ]  +  V  +  Ik/U  +  kn),Juv}  +  Pg/Cl  +  kn) 


=  [e^/(l  +  kn)  ^  Cr  +  n  coS(())'^]Kl  +  kn)  ^  (r  +  n  coS(|))'^t] 

(2) 


n 


where 


T  = 


y [u^  -  ku/ (1  +  kn) ] 


(3) 


n-Momentum  Equation; 

p{u[Vg/(l  +  kn)]+  -[k/(l  +  kn)]u^}  +  ^  C4) 

where  with  the  thin  shock  layer  approximation,  equation  C4)  becomes, 

2 


^n  ^  u 


(5) 


Energy  Equation; 


p{u[Tg/(l  +  kn)3+  V  Tj^}  -  {u[p^/(l  +  kn)  ]  +  v  p^^} 

=  [e^/(l  +  kn)  (r  +  ncos(j))'^]  [  (1  +  kn)  (r  +  ncosiji)  “^q] . 

,  /  -2  ,  .  2 
+  (e  /y)T 


(6) 


where  q  =  (y/a)  Tj^  (.7) 

Equation  of  State; 

P  =  [(y-1)/y]Tp  (8) 

Viscosity  Law 

y  =  [(1  +  c')/(T  +  c') (9) 
c'  =  c*/(y  -  DMf  T*  (10) 


* 

where  c  is  taken  to  be  198. 6°R  for  air. 
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Surface  slip  conditions  consistant  with 
the  above  set  of  equations  are  given  as 

V  =  0 

u  =  a3^(l/p){  [(y-1)/y]T}^/2t 

p  =  (a/T) { [ (y-1)/y]T}^^^ 

T  =  C3^(a/p)  {  [  (y-1)/y]T}^/^ 


the  approximations  used  in 


(11a) 


k  at  n  =  0 


(11b) 


q 


(11c) 


q  J 


(lid) 


where  t  and  q  are  given  by  equations  (3)  and  ,C7)  . 

Conditions  at  the  shock  surface  are  obtained  using 
the  concept  of  "shock  slip"  to  represent  the  usual  higher  order 
Reynolds  number  effects  on  the  shock  compression  process.  This 
gives  modified  Rankine-Hugoniot  relations  as  shown: 


Ugh  =  sin(a+e)  +  v^j^  cos(a+e) 

Vgj^  =  -  Ugj^  cos(a+e)  +  Vgj^  sin(a+B) 


(12) 

(13) 


1  1 

where  u^j^  and  v^j^  are  the  components  of  velocity  tangent  and 

normal  to  the  shock  interface,  respectively,  and  are  given  along 
with  the  temperature,  pressure  and  density  from  the  following 
expressions : 

Psh  ^Ih  =  - 


2 

£ 


U  .  (ul)  ,  + 
•^sh  n  sh 


sino.  Ugj^ 


sina  cosa 


(15) 


Psh  sh  ^  ^sh  -(sina/2)  Cu^^^  -  cosa)  ^ 

=  sina/2  {  [4y/(y+1)  ^Isin^a  +[2/(y-11  -  4  (y-1)  /  (  y+D  ^  ]  l/M^ 
-  4/(y+1)^  wf  sin^a)  (16) 
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Pgj^  =  [2/(y+1)  ]sin^a  -  (y-l) /y  (y+DM^ 
•^sh  ^ 


(17) 

(18) 


For  reasons  discussed  in  Reference  ( l  )  the  above  equations 
were  solved,  here  after  the  independent  and  dependent  variables 
were  normalized  by  their  corresponding  shock  values.  The 
resulting  set  of  governing  equations  and  the  boundary  conditions 
are  given  in  Reference  (  1  ) . 


III.  NUMERICAL  ANALYSIS 

This  set  of  viscous  shock  layer  equations  is  of  the 
parabolic-hyperbolic  type  and  is  therefore  solved  using  a 
method  similar  to  that  used  for  solving  the  boundary  layer 
equations,  such  as  the  method  of  Blottner  and  Flugge-Lotz 
(Ref.  10).  The  calculation  of  the  shock  shape,  though, 
represents  an  elliptic  effect.  To  properly  account  for  the  shock 
shape,  the  present  method  of  solution  adopted  a  relaxation 
technique  wherein  the  shock  slope  was  determined  iteratively. 

After  each  pass  along  the  body  surface  a  least  squares 
Chebyshev  polynomial  was  used  to  smooth  the  numerically 
calculated  shock  shape  over  the  entire  range  of  integration 
in  the  downstream  direction.  This  method  of  solution  was 
found  to  work  well  for  flow  past  blunt  bodies  where  the  shock 
stand  off  distance  is  not  small  and  the  shock  shape  differs 
significantly  from  the  body  shape.  For  a  detailed  description 
of  the  numerical  method  one  is  referred  to  Reference  ( 9  ) . 

The  overall  method  of  solution  employed  is  as  follows ; 

An  initial  guess  was  made  on  the  shock  slope.  This  was 
done  by  giving  initial  values  to  the  coefficients  of  the 
Chebyshev  polynomial  which  represents  the  shock  slope  function. 

The  equations  were  then  solved  by  starting  with  the  stagnation 
point  where  they  reduce  to  a  set  of  ordinary  differential 
equations.  The  first  equation  solved  was  the  energy  equation 
so  that  thereafter  all  quantities  such  as  viscosity  related  to 
temperature  could  be  evaluated^  Next  the  s-momentum  equation 
was  integrated  to  determine  a  u  velocity  profile,  and  then  the 
continuity  equation  was  solved  to  determine  first  the  shock  stand 
off  distance  from  equation  (39a)  (Ref.  9  )  ,  and  then  the  v  component 
of  velocity  from  equation  (31)  (Ref . 9  ),  Finally  equation  C33) 

(Ref.  9  )  was  integrated  to  determine  the  local  pressure  level. 
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Repetition  of  the  above  steps  at  a  given  station  continued 
until  the  solution  converged-  The  method  then  stepped  along 
the  body  surface  and  iterated  at  each  station  to  achieve 
converged  solutions.  To  accelerate  the  convergence  process^ 
the  previous  station  values  of  the  profiles  were  used  at  each 
new  step  as  a  first  guess.  Once  the  above  method  had  passed 
over  the  entire  mesh,  the  coefficients  of  the  Chebyshev 
polynomials  were  recalculated  using  the  calculated  values  of  the 
shock  derivative  at  every  station..  The  entire  procedure  was 
repeated  until  the  coefficients  of  the  Chebyshev  polynomial 
converged  to  a  desired  accuracy. 

The  solutions  so  obtained  were  labeled  as  ‘^Thin  Shock 
Layer"  solutions  when  the  approximate  normal  momentum 
equation  (Eqn.  5)  was  used.  Further  solutions  using  the  more 
correct  normal  momentum  equation  (Eqn.  4)  were  obtained  by  a 
similar  procedure.  These  solutions  were  labeled  as  "Full 
Shock  Layer"  solutions.  However  for  this  phase  of  calculations 
the  V  terms  in  the  normal  momentum  equation  were  used  from  the 
previous  iteration . 


The  shock  slip  conditions  were  handled  by  evaluating  the 
(u^)gh  ^"^n^sh  from  the  previous  step  in  the  Iteration 

and  then  solving  the  resulting  equation  for  this  shock  condition. 


IV.  RESULTS  AND  DISCUSSION 

Figure  (2)  presents  a  comparison  of  normalized  stagnation 
point  heating  on  spheres  with  experimental  data  and  other 
theoretical  results.  The  numerical  calculations  were  made  at  a 
free  stream  Mach  number,  of  8,  a  wall  to  stagnation 

temperature  ratio  of  0.25  and  a  total  temperature  of  500 ®R. 

The  heat  flux  in  figure  (2)  is  normalized  with  its  value 
obtained  from  a  boundary  layer  analysis  at  the  same  test 
conditions.  Note  that  the  full  layer  results  rightfully 
approach  unity  as  Reynolds  number  increases  in  good  agreement 
with  the  higher  order  boundary  layer  theory  of  Van  Dyke 
(Ref.  21)  that  gives 


q/qg  1  +  0. 866//Reg 


Comparison  of  the  present  full  shock  layer  results  Cincluding 
shock  and  wall  slip  effects)  with  the  experimental  data  of 
Potter  and  Miller  (Ref.  11),  Wittliff  and  Wilson  (Ref.  12) 
and  Hickman  (Ref.  13)  is  good  over  the  entire  Reynolds  number 
range . 
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Figure  (2)  also  shows  the  analytical  result  of  Cheng 
(Ref»  6  )  where  a  Newtonian  thin  shock  layer  approach  was 
employed*  Both  wall  and  shock  slip  effects  were  incorporated 
in  these  solutions  and  thus^  are  seen  to  compare  favorably 
with  the  present  calculations  over  the  entire  Reynolds 
number  range  * 

In  an  effort  to  assess  the  importance  of  wall  and  shock 
slip  conditiohs  on  the  predicted  results^  solutions  were 
also  obtained  without  slip  effects  included*  As  shown  in 
figure  (2)  as  the  Reynolds  number  decreases ^  wall  slip 
effects  become  important  and  then  shock  slip  effects  grow 
to  a  significant  level*  Both  these  effects  are  of  equal 
significance  in  the  low  Reynolds  nurober  regime* 

It  is  also  noted  from  figure  (2)  that  while  the  full 
shock  layer  equations  yield  results  that  compare  well  with 
various  data^  their  thin  layer  approximations  give  results 
that  are  significantly  lower  than  the  data  for  the  Reynolds 
number  range  shown  * 

However^  note  should  be  made  that  the  data  of  Ferri 
et*  al*  (Ref.  14)  which  were  obtained  with  M  -  8*  -  0*25 

and  “  2300®R^  was  found  to  be  higher  than  those  predicted 

by  any  of  the  other  analytical  or  experimental  results* 

Hickman  (Ref*  13)  attributed  this  apparent  increase  to  the 
higher  stagnation  temperature  level  of  Ferri’s  experimental 
conditions  (2300° R)  as  compared  to  his  test  stagnation 
temperature  level  of  500°R*  In  order  to  clarify  this  situation 
another  set  of  numerical  calculations  were  made  where  the  test 
conditions  were  taken  to  be  the  same  as  those  used  by  Ferri 
et*  al*  (Ref*  14)*  The  results  of  such  a  calculation  is  shown 
in  Figure  (3)*  It  is  noted  from  this  figure  that  whan  the 
stagnation  temperature  of  the  test  condition  is  increased  from 
540°R^  as  used  by  Hickman  (Ref*  13)^  to  2300°R  corresponding 
to  the  test  case  of  Ferri  et*  al*  (Ref*  14) ,  the  stagnation 
point  heating  levels  do  not  show  much  influence*.  Hence  the 
difference  between  the  data  of  Ferri  et*  al*  (Ref*  14)  and  the 
present  calculations  remains^  as  yet^  unexplained* 

In  order  to  further  test  the  viscous  shock  layer  model , 
numerical  calculations  were  made  corresponding  to  the  theoretical 
calculations  of  Dellinger  (!Ref*,  5  ).  at  a  free,  stream  Mach- 

numfoer,  of  14*3,  a  free  stream  Reynolds  number,  Re  ,  of 

3200  and  a  wall  to  stagnation  temperature  ratio  of  0*o282* 

Figure  (4)  presents  the  shock  wave  and  sonic  line  locations  for 
this  case*  The  results  fx"om  the  present  analysis  compare  well  with 
the  calculations  of  Dellinger  and  with  the  inviscid  blunt  body 
solution  of  Van  Dyke  and  Gordon  (Ref*  15)*  The  apparent 
descrepancy  in  the  location  of  the  sonic  line  between  the  present 
results  and  the  inviscid  results  of  Van  Dyke  and  Gordon  (Ref*  15) 
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near  the  shock  is  due  to  the  fact  that  the  present  results  are 
for  a  Mach  number  of  14.3  as  against  10  used  by  Van  Dyke  and 
Gordon  (Ref.  15 ) •  The  difference  near  the  body  is  due  to  the 
viscous  effects  which  cause  the  Mach  number  to  be  lower  than 
inviscid  theory  would  predict.  Thus  the  results  obtained 
at  this  test  condition  are  seen  to  be  reasonably  good. 

Further  test  data,  such  as  those  due  to  Boylan  CBef.  16 )  are 
available  at  the  higher  free  stream  Mach  number  of  21.0,  than 
those  shown  in  Figure  (2)  .  Figure  (5)  shows  the  Stanton 
number  at  the  stagnation  point  corresponding  to  the  test  data 
(M^  =  21,  =  0.11)  of  Boylan  (Ref.  16),  The  theoretical 

calculations  shown  in  this  figure  were  obtained  using  the  "thin 
layer"  version  of  the  full  shock  layer  equations.  It  is  noted 
that  these  equations  correctly  predict  the  qualitative  behavior 
of  the  stagnation  point  Stanton  number  with  increasing  altitude, 
however  the  quantitative  values  are  seen  to  be  significantly 
different  from  those  of  the  data.  This  difference  is  seen 
in  Figure  (5)  to  reduce,  first  with  the  inclusion  of  the  effects 
of  wall  slip, and  then  further  with  the  inclusion  of  the  effect 
of  shock  slip.  The  influence  of  the  slip  is  seen  to  be  smaller 
at  lower  altitudes,  as  seen  also  in  Figure  (2),  Hence  the  "thin 
layer"  version  of  the  full  shock  layer  equations  would  not 
predict  correct  physical  properties  at  lower  Reynolds  numbers. 
However,  as  seen  in  Figure  (6),  the  full  layer  equations  without 
slip  effects  included  are  observed  to  give  as  erroneous  results 
as  their  "thin  layer"  counterpart.  These  equations  are  seen 
to  give  reasonably  good  results  only  when  both  wall  and  shock 
slip  effects  are  included.  Thus,  the  effects  of  slip  are  seen 
to  be  large  in  the  lower  Reynolds  number  regime  and  must  be 
accounted  for. 

In  order  to  assess  the  capability  of  the  present  model 
to  represent  the  full  flow  structure  acrosrs  the  ihock  layer, 
comparison  is  made  in  Figure  (7)  with  the  density  profile  measured 
by  Russel  (Ref.  17)  for  flow  of  nitrogen  past  a  sphere  at 
various  free  stream  Reynolds  numbers.  Figure  (7a)  shows  the 
density  profile  corresponding  to  a  shock  Reynolds  number  of  129 
and  a  free  stream  Mach  number  of  4.11.  The  experimental  result 
shows  a  finite  thickness  of  the  shock  wave  whose  effect  on  the 
shock  layer  density  is  seen  to  be  correctly  modeled  using  the 
present  approach.  It  is  noted  that  the  shock  discontinuity 
modeled  in  the  present  analysis  rightly  occurs  in  the  middle 
of  the  shock  thickness  shown  by  the  experimental  data..  This 
figure  also  shows  the  results  of  Jain  and  Adimurthy  (Ref.  18) 
where  solutions  were  obtained  using  a  series  truncation  technique 
of  the  Navier  Stokes  equations  in  the  stagnation  region..  While 
their  results  well  model  the  shock  layer  region,  they  apparently 
misrepresent  the  shock  wave  region  due  to  the  fact  that  the 
series  truncation  approach  tacitly  assumes  the  shock  to  be 
spherical  over  the  entire  stagnation  region.  The  "thin  layer" 
version  of  these  equations  are  also  seen  to  yield  results  that 
compare  well  with  the  data«.  It  is  noted  from  this:  figure  that  the 
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effects  of  wall  and  shock  slip  are  not  large  so  far  as  comparison 
with  density  profile  data  is  concerned^  however  these  effects 
are  seen  to  cause  a  significant  change  in  the  wall  properties, 
Thus^  care  must  be  taken  to  include  these  effects  in  the  proper 
regime. 

Figure  (7b)  shows  the  density  profile  corresponding  to 
Russel's  data  of  a  free  stream  Mach  number  of  4ol9  and  a  shock 
Reynolds  number  of  50,  These  results  show  that  with  the 
decreasing  shock  Reynolds  number  the  density  profile  predicted 
by  "thin"  and  "full"  shock  layer  equations  are  significantly 
different.  Further  decrease  in  shock  Reynolds  number  causes 
a  large  difference  between  the  density  profile  predicted  by  the 
full  shock  layer  equations  and  their  "^thin  layer counterpart 
as  shown  in  Figure  (7c)  where  the  test  conditions  were  taken  to 
be  a  free  stream  Mach  number,  of  4,38  and  a  shock  Reynolds 
number  of  24,0,  The  experimental  data  in  this  figure  shows  a 
comparatively  larger  shock  wave  region  than  that  seen  in 
Figures  (7a)  or  (7b) ,  however  the  viscous  shock  layer  model 
predicts  the  density  profile  reasonably  well.  The  slip  effects 
are  noted  to  have  a  significant  influence  on  the  predicted 
density  profile  at  this  low  shock  Reynolds  number  case.  Thus 
it  is  observed  from  these  figures  that  when  the  full  shock  layer 
equations  are  used  including  the  effects  of  slip,  the  predicted 
density  prof ile  compares  well  with  the  experimental  data  for 
a  wide  range  of  shock  Reynolds  numbers » 

The  density  profiles  discussed  above  were  obtained  for  a 
diatomic  gas  such  as  nitrogen  for  which  the  ratio  of 
specific  heats  were  taken  to  be  1,4,  Further  test  data  were 
also  considered  for  a  monotomic  gas  such  as  argon  where  the 
ratio  of  specific  heats  were  taken  to  be  1,667,  Figures  (8a,b,c) 
show  the  density  profiles  for  such  a  gas  for  various  free  stream 
Reynolds  numbers  corresponding  to  the  test  condition  of 
Russel  (Ref,  17),  Figure  (8a)  shows  the  density  profile  for  a 
shock  Reynolds  number  of  358  and  a  Mach  number  of  3,94,  The 
predicted  results  compare  well  with  the  data  and  the  effects 
of  thin  layer  approximations  and  slip  are  not  significant  in 
a  region  away  from  the  wall.  However  these  effects  grow  to  a 
significant  level  as  the  shock  Reynolds  number  is  decreased 
to  a  value  of  96,0,  as  shown  in  Figure  (8b) ^  and  further  to  a 
value  of  37,0  as  shown  in  figure  (8cl,  No  different  outcome 
thus^  is  observed  when  the  gas  model  is  changed  from  a  diatomic 
gas  to  a  monotomic  gas.  As  a  final  comparison,  Figure  (9)  shows 
the  temperature  profile  on  the  stagnation  line  of  a  spherical 
body  at  M  =  19,  T  /T  =0.19  and  Re  =  54,0  measured  by 

Ahouse  and  Bogdonoff  (Ref,  19)  along  with  various  other 
theoretical  calculations.  The  results  obtained  by  the  present 
analysis  only  compare  qualitatively  with  the  experimental  data, 
but  they  do  show  a  trend  similar  to  the  shock  layer  analysis  of 
Cheng  (Ref,  6  )  and  the  Navier““Stokes  analysis  of  Li  (Ref,  20)  , 


12 


AEDC-TR-76-53 


Note  that,  even  though  the  similarity  solution  of  Jain  and 
Adimurthy  (Ref.  18)  tacitly  assumes  the  shock  to  be  spherical 
in  the  nose  region,  their  results  seem  to  be  surprisingly  closer 
to  the  experimental  data.  The  exact  source  of  the  discrepancy 
between  the  present  shock  layer  analysis  and  the  experimental 
data  in  the  outer  reaches  of  the  layer  £s:  as^  yet  unexplained.. 

In  an  overall  sense,  the  results  obtained  here 
indicate  that  the  viscous  shock  layer  model  compares  well  with 
experimental  data  to  very  low  Reynolds  numbers.  It  is  found 
that  with  the  inclusion  of  the  shock  and  body  slip,  the  full 
viscous  shock  layer  model  apparently  enjoys  a  range  of  validity 
down  to  shock  Reynolds  numbers  on  the  order  of  20-30. 


V.  CONCLUSIONS 


A  finite  difference  method,  where  solutions  are  obtained 
using  a  least  square  Chebyshev  polynomial  to  fit  the  numerically 
determined  shock  shape  (Ref.  7) ,  has  been  successfully  applied 
to  the  solution  of  the  viscous  shock  layer  equations  for  flow 
past  a  sphere.  Solutions  so  obtained  were  used  to  establish 
the  range  of  validity  of  these  equations  by  comparison  with 
various  experimental  data.  Based  on  these  comparisons  it  is 
believed  that  the  full  shock  layer  model  when  used  along  with 
the  wall  and  shock  slip  conditions  apparently  enjoy  a  range  of 
validity  down  to  shock  Reynolds  numbers  on  the  order  of  20  to 
30.  The  calculated  surface  heating  levels  and  the  shock  layer 
profiles  are  seen  to  be  in  good  agreement  with  the  various 
experimental  and  theoretical  results.  The  approach 
in  the  present  analysis  is  seen  to  well  represent  the  shock 
wave  location  and  other  physical  quantities  of  the  flow.  It  is 
further  determined  that  the  thin  layer  version  of  these  equations 
are  inadequate  for  low  values  of  the  shock  Reynolds  numbers. 
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Figure  1.  Coordinate  system. 
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Figure  3.  Further  comparison  of  stagnation  point  heat  transfer. 
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Figure  4,  Shock  wave  and  sonic  line  locations  for 
a  spherical  body. 
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Figure  5.  Variation  of  Stanton  number  with 
Reynolds  number  for  thin  layer. 
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Figure  9 .  Stagnation  point  temperature  profile • 
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LIST  OF  SYMBOLS 


Body  nose  radius  of  curvature 

Slip  constant  taken  to  be  1. 2304 (2-6^) /0^ 

Slip  constant  taken  to  be  1. 1750  (2-e^)/e^ 

Slip  constant  taken  to  be  2 . 3071 (2-a^) /a^ 

Surface  curvature 
Free  stream  Mach  number 

Coordinate  measured  normal  to  the  body,  nondimensionalized 
by  the  body  nose  radius 

*  2 
Pressure  p  /(p  U  ) 

Nondimensional  axisymmetric  radius  measured  to  a  point 
on  the  body  surface 

*  it  it  it 

Shock  Reynolds  ntimber,  a 

Defined  as  Re^/t/2  ,  e  =  0*13 

it  it  it  it 

Defined  as  p  U  a  /y^ 

^00  00  O 

Surface  distance  coordinate  measured  along  the  body 

*  ,  ,  itO  *4 

Temperature,  T  =  T  / (U  /C  ) 

P 

Free  stream  temperature 

*  it 

Velocity  component  tangent  to  the  body  surface,  Cu  /U^) 
Free  stream  velocity 

Component  of  velocity  tangent  to  shock  interface 

iSf  it 

Velocity  component  normal  to  the  body  surface,  (v  /U^) 

Component  of  velocity  normal  to  shock  interface 
Shock  angle,  see  Figure  1 

Thermal  accommodation  coefficient  here  taken  to  be  1 
Angle  defined  in  Figure  1 
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Symbol 


Ratio  of  specific  heats 

a  *2  *v  *  *  *  1/2 

Perturbation  parameter,  e  =  [y  (U^.  /C  )/p^  a  ]  ' 
Fraction  of  incident  molecules  diffusely  reflected 

i<  *  *2  * 

Coefficient  of  viscosity,  y  =  y  /y  (U  /C) 

00  p 

* 

Density,  p  =  p  /p^ 

Free  stream  density 

*  *  *2 
Shear  stress,  t  /(p„  ) 

Body  angle  defined  in  Figure  1 

Prandtl  number 


Subscripts 

1  Wall  value 

0  Stagnation  point  value 

sh  Behind  the  shock 

“  Free  stream  conditions 


Superscripts 

*  Dimensional  quantities 

J  0  for  plane  flow  and  1  for  oxlsymmetric  flow 

(”)  Dependent  variables  normalized  with  their  respective 

shock  values 
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